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Abstract. Surface diffusion of interacting adsorbates is here analyzed within 
the context of two fundamental phenomena of quantum dynamics, namely the 
quantum Zeno effect and the anti-Zeno effect. The physical implications of 
these effects are introduced here in a rather simple and general manner within 
the framework of non-selective measurements and for two (surface) temperature 
regimes: high and very low (including zero temperature). The quantum 
intermediate scattering function describing the adsorbate diffusion process is then 
evaluated for flat surfaces, since it is fully analytical in this case. Finally, a 
generalization to corrugated surfaces is also discussed. In this regard, it is found 
that, considering a Markovian framework and high surface temperatures, the anti- 
Zeno effect has already been observed, though not recognized as such. 



1. Introduction 

Traditionally, one finds in the literature that activated surface diffusion and related 
phenomena, such as the frustrated translational motion (the so-called T-mode), are 
well described by classical models whenever heavy particles and long time scales (of 
the order of tens and hundreds of picoseconds) are considered. In such cases, which 
are the typical ones, any trace of quantumness is essentially swept out; one only finds 
clues of the lost quantumness through "reminders" , for example, in recoil energies. 
However, below these time scales also explored through collisions between the probe 
particle and the diffusive one, quantum effects should be stronger. In this domain, the 
commutation rule for position operators at different times should play an important 
role in the corresponding dynamics, for example. In principle, this time domain would 
be ruled by a coherence time tc ~ ^/3(/3 = l/ksT and ks is the Boltzmann constant), 
determined essentially by the surface (substrate) temperature T. Thus, in general, 
for times smaller than tc-, quantum effects should be dominant over thermal effects. 
Hence, a question that arises in a natural way is how and why this transition takes 
place or, in other words, why the diffusing, interacting adsorbates lose their quantum 
behavior. Within the scenario depicted by this quantum-to-classical transition, one 
can understand surface diffusion processes as decaying processes, with their evolution 
being monitored along time through the so-called intermediate scattering function. 
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Transitions like the aforementioned one are strongly connected with measure 
processes in the sense that, regardless of how a measure is defined, in general such 
processes can be understood as a system coupled to an external agent acting on it 
continuously or at certain times. With each interaction between the system and the 
measuring device, the former would lose part of its quantumness, thus leading to 
typical classical behaviors (decoherence) . This brings us to the scenario of the quantum 
Zeno effect (QZE) and its complementary, the anti-Zeno effect (AZE), which consist 
of precisely inhibiting or enhancing, respectively, the system natural (quantum) decay 
by carrying out a series of measures on it. These effects appear and are analyzed 
commonly within the context of unstable quantum systems. However, here we would 
like to introduce them within a very different context, namely surface diffusion of 
interacting adsorbates. In this regard, a novel analysis and (to some extent) new 
interpretation of quantum surface diffusion in terms of QZE and AZE (excluding 
tunneling mediated diffusion) is presented. The primary motivation of this work 
arises from the possibility at present to carry out surface diffusion experiments finely 
resolved in time by means of neutron [1] and '^He spin-echo techniques [2j|3]. From 
these experiments, one obtains the time-evolution of a quantity (namely the so-called 
polarization function) proportional to the intermediate scattering function, 

I{AK,t) = (e-»AK.R(0)g^AK.R(^)^_ 

This time-dependent (correlation) function is the space Fourier transform of the time- 
dependent pair correlation function or G-van Hove function, G(R, i), a generalization 
of the pair-distribution function from the theory of liquids ■ Given an adparticle is 
at some position R(0) on the surface at time t = 0, G(R, t) accounts for the averaged 
probability of finding the same or another adparticle at R at a time t (here, R denotes 
the adparticle position operators given in the Heisenberg picture). To some extent, 
one could therefore interpret this function as a time decay law, describing the loss 
of correlation between a system initial configuration, represented by e*^^'^(°\ and 
its configuration at time t, e*^^'^(*\ measured along the AK-direction (AK is the 
wave- vector transfer parallel to the surface of the probe particles, i.e., neutrons or ^He 
atoms) . 

Experimentally, it is observed that /(AK,t) decays smoothly with time, 
quadratically at short times and like e~ at longer times, with rjcS depending on 
the reciprocal lattice associated with the surface, the momentum transfer and the 
friction coefficient r/ ruling the diffusion time scale. In activated surface diffusion, 
the quantum Langevin equation formalism has been considered to evaluate the 
intermediate scattering function ([ij. The Caldeira-Leggett Hamiltonian model has 
also been used to describe the diffusion dynamics, representing the surface as an 
infinite collection of harmonic oscillators at a given surface temperature T (reservoir) 
and with an Ohmic (constant) friction (one-bath model [B]). If interacting adsorbates 
are considered, the total friction splits up into two contributions, ?7 = 7 -I- A, as shown 
by means of the so-called interacting single adsorbate (ISA) approximation, which has 
been applied with success at low and intermediate coverages (two-bath model 7, 8J. 
Within this model, 7 is associated with the substrate or reservoir, while A represents a 
collisional friction due to the collisions among adsorbates, thus being connected with 
the surface coverage. By invoking the elementary kinetic theory of transport of gases, 
a simple relation can be found between the collisional friction and the coverage 9 at a 
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temperature T 




6p9 keT 



(2) 



where a is the unit ceU length of an assumed square surface lattice and p is the effective 
radius of an adparticle of mass m. From a quantum- mechanical viewpoint, when the 
adsorbate interacts with the surface the entanglement with the degrees of freedom of 
the latter leads to the loss of coherence (or decoherence) of the former, i.e., its ability to 
display quantum behaviors in the corresponding dynamics (e.g., quantum interference) 
[5] . In the context of quantum Brownian motion [TU] , this phenomenon where a highly 
delocalized state in position and/or momentum transforms into a localized classical 
state is called environment induced decoherence. In the particular case of the adsorbate 
surface diffusion, the role of the measurement process on the adsorbates would be 
associated with the surface and the other surrounding adsorbates. This process thus 
belongs to the class of non-selective, continuous and indirect measurements, which 
are different to the more common projective (von Neumann) measurements. As a 
consequence, the adsorbate position would be describable by a Langevin equation, 
where the apparatus would be included through an interaction Hamiltonian llj . More 
specifically, the apparatus would be a reservoir [T21[T3], where the effect of indirect, 
non-selective measurements (e.g., collisions) would contribute to the dephasing of the 
system quantum state, making the phase of the adparticle completely random [14]. 
Within the two-bath model, the surrounding adsorbates form a reservoir, such that 
their collisions with the tagged adsorbate (and therefore their frequency, measured 
through A) can be controlled through the coverage after Eq. 

Finally, we would like to mention that in the context of atom-surface scattering, 
Levi has recently introduced and discussed decoherence [H] and the QZE [16 . 

The organization of this work is as follows. To be self-contained, in Section [2] the 
physics behind the QZE and AZE is introduced in a brief manner. In Section [3] an 
overview of the wave packet diffusion dynamics through flat surfaces is presented in 
order to introduce some concepts and notions on QZE and AZE within the context of 
surface diffusion. We would like to stress here that the appearance of QZE and AZE 
is very much related to the type of measurement carried out, leading to conclusions 
which can be quite different. In Section |4l the intermediate scattering function for 
flat and corrugated surfaces is analyzed, seeing it as a decay rate in the Heisenberg 
picture, and provides a new interpretation of the experimental and theoretical results 
previously reported. Finally, in Section [SJ some conclusions are present together with 
possible extensions of this type of analysis to the frustrated translational mode. 

2. Elementary grounds of QZE and AZE 

In the standard fashion, the QZE and AZE are usually introduced within the context 
of a series of ideal, repeated measurements (or projections) on a system during its time- 
evolution [T3 l [T4 l[T7H26] . Thus, let us first focus on the short-time decay dynamics of 
an unstable system described by the total Hamiltonian 



where accounts for the system free dynamics and V describes some interaction 
that ultimately leads to the system decay. At t = 0, the system is supposed to be 
in a pure state \'lpo)^ which is a normalized vector in Hilbert's space. The survival 



H = Ho + V, 



(3) 
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probability or probability to find the system in the same state at a certain subsequent 
time t is given by 

P{t) = \{^Jo\m)\^ = |(V'oIe-^^*/^|Vo)r (4) 

At very short times, P{t) can be expressed by Taylor expansion of the evolution 
operator up to second order as 

F(t)«l-((^o|^r2|V.o)-(V'oIi?|V'o)');^ = l- . (5) 

That is, at very short times the survival probability decays quadratically with time. 
If initially the system is in an eigenstate of Hq , then we find the characteristic Zeno 
time Tz depends essentially on the variance of the interaction potential, i.e., 

^ ^ /g^, 

V{MV'\r^o)-{MV\^o)^' 

Now, consider a number N of instantaneous ideal measurements (projections) are 
performed at very small intervals of time r, such that t = Nt, in order to ascertain 
whether the system still remains in its initial state or not. Each time a measurement 
is performed, the system wave function is "collapsed" and its subsequent evolution 
starts again from the state |'0o)- Therefore, the probability to find the system in this 
state after a time t and after performing N measurements will be 

OH N 



PW(t)=[P(T)] 



N 



1 ( ' 



(7) 



\Ntz 

In the limit N ^ oo, this quantity becomes 

« e-*'/^^z _^ 1. (8) 

That is, in the ideal limit where the system can be monitored indefinitely, its quantum 
state will remain the same, without evolving, because the interaction is not enough to 
remove its quantum coherence. In this way, by means of consecutive measurements, 
the system decay is slowed down. This is the QZE, predicted by Misra and Sudarshan 
in 1977 for an unstable particle [T7]. Later on, in 1989, Itano et al. [32] observed how 
Be atoms were inhibited to evolve into the excited state by means of ultraviolet pulses. 
More recently, Fischer et al. [24] observed QZE with cold Na atoms in an experiment 
trying to mimic the proposal of Misra and Sudarshan. 

In the experiments conducted by Fischer et al. [24], AZE was also found, which 
was predicted by Kofman and Kurizki |141I23] within the density matrix formalism. 
This phenomenon occurs when the decay process is accelerated due to continuous 
measurements (ideal and sufficiently frequent). Notice that the survival probability 
([7]) can also be expressed as a general function of the measurement time, r, as [25ll26] 

F(^)(t) = e^''""-^^)*, (9) 

If r is small compared to tz (although not always this constitutes an appropriate time 
scale [26]), then ~^ t/t"^, as in ([8]). However, for t large enough (but still small), 
the system should display a typical decay at a constant rate 7froc, as it is found in 
unstable systems according to Fermi's golden rule, i.e., 7off — > 7fioc- The accelerated 
decay with respect to the exponential decaying behavior is the AZE. The general decay 
law ruling the behavior of unstable quantum systems has been discussed in detail by 
Peres [27]. The physical conditions for leading to the observation of QZE and AZE 



Quantum Zeno and anti-Zeno effects in surface diffusion 



5 



have been established in a very illuminating way by Kofman and Kurizki |14| . If the 
spectral density of states coupled to the initial state is a dense band or continuum 
(acting as a reservoir), the measurement-modified decay rate can be expressed as 
a simple overlap between the spectral density of final states or reservoir coupling 
spectrum, G'(a;), and the measurement-induced initial state level width, F{uj\t), 
according to 



where the measurement frequency A = r^^ (we use the bar over A in order to 
avoid confusions with the coUisional frequency from surface diffusion), is related to 
the initial state energy uncertainty, Ai?, as AE/X ^ h. Hence the decay rate is 
essentially determined by the spectral density profile within a bandwidth around its 
energy level. According to the same authors [14], pO)) constitutes a universal result: 
frequent measurements on a given initial state generally lead to its dephasing through 
randomization of the corresponding phase. The broadening of the initial state is of the 
order of the measurement frequency and it can be seen as an analog of the coUisional 
broadening leading to a phase randomization of the state. Two extreme cases are then 
envisaged. On the one hand, when the measurement frequency A is much greater than 
the spectral density width and the detuning between the reservoir center of gravity and 
the initial state frequency position, the QZE holds (the decay rate goes like A^^). A 
reduction of the decay rate is then obtained when compared to the measurement-free 
(or Fermi's golden rule) decay. Mathematically, the spectral density is assumed to be 
a Dirac S peak. In other words, when an Ohmic friction or white noise, characterized 
by a linear spectral density, is assumed, the QZE will not be observable. Non-Ohmic 
reservoirs should be then considered with a cut-off frequency. On the other hand, if the 
measurement frequency is much smaller than the corresponding detuning, the decay 
rate is shown to grow with A, which leads to AZE. Actually, the AZE seems to be 
much more ubiquitous than the QZE [14 . More recently, Maniscalco et al [55], within 
the density matrix formalism, have established the conditions for the occurrence of 
such effects within the quantum Brownian motion; in particular, they have studied 
a quantum harmonic oscillator linearly coupled to a quantum reservoir modeled as 
a collection of non-interacting harmonic oscillators at thermal equilibrium. During 
the time-evolution, the system was subject to a series of non-selective measurements, 
i.e., measurements which do not select the different outcomes 11 . The factorization 
displayed in ([7]) follows from the fact that at second order in the coupling, the density 
matrices of the system and environment factorize at any time (25] . 

Within the scenario of the surface diffusion of interacting adsorbates, and 
according to PH)) . the coUisional friction A would govern the F function, and the 
surface friction 7, the G function. Furthermore, it is interesting to note the following 
comparison. Within the standard projective-measurement scenario described above, 
a series of N measurements is carried out at regular intervals of time; within the non- 
selective-measurement scenario, A provides an average number of collisions per time 
unit [30j and, therefore, an average number of measures per time unit, although such 
measures are carried out at random (according to a Poissonian distribution). 




(10) 
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3. Surface diffusion of a wave packet by a flat surface 

Consider the diffusion of a wave packet through a fiat surface. It is the simplest 
integrable dissipative quantum system we may devise. This simple example is quite 
illustrative to show the different conclusions about the appearance of the QZE and 
AZE when carrying out ideal (or projections) and indirect measurements. 



3.1. High surface temperatures 

Within the ISA model, the quantum Langevin equations for the Heisenberg position 
operators read (for Ohmic friction or linear spectral density) as 

x(t)^-T^x{t)+5F,{t), 

m = -vy{t) + SFyit), ^''> 

where R = (x, y) denotes the position operators of a single adsorbate moving on a 
flat surface and the "dots" over the position operators denote time-derivatives. In 
this model, two non-correlated noises (per mass unit) are assumed to simulate the two 
baths: a Gaussian white noise, accounting for the lattice vibrational effects due to the 
surface temperature and leading to the interaction with the adsorbates, and a white 
shot noise, which simulates the adsorbate-adsorbate collisions. Then, for each degree 
of freedom, we have 6Fi{t) = 5Fi^a{i) + ^Fi,s{t) where i — x,y, the noise fluctuation 
is given hy SF = F — (F), and G and S stand for the Gaussian and shot noise, 
respectively. At high surface temperatures, /3~^ » Hrj (or r]~^ ^ tc ^ hfi), noise is 
considered to be classical and its autocorrelation function at two different times is well 
described by a Dirac (5- function, as assumed in ()lip . In other words, the time interval 
between these two times has to be greater than tc. 

Consider an adparticle of mass m initially placed on a given position of a flat 
surface and represented by a Gaussian state 

V'(x, y, 0) - -L= e-^^'^'^o-vV^'^l , (12) 

where the initial width along each direction is the same and equal to ctq- The adparticle 
is assumed to be initially in equilibrium with the reservoir or heat bath (surface) at 
a temperature T, but weakly coupled to the environment that dissipation can be 
neglected. The role of the initial conditions has been very often discussed in the 
literature (see, for example, [EUSI])- After a time t, the probability to find the 
particle at a given position (x, y) is given by averaging the survival probability over a 
thermal (Maxwell-Boltzmann) distribution of velocities [12,31-35j 

P{x,y^t) = ^——^ e-V2»^W-yV2«.?(*). (13) 

2TTWx{t)Wy{t) 

Note here that the interaction with the environment makes the quantum state 
describing the system to pass from pure to a statistical mixture. Therefore, the 
probability ([T^ has to be interpreted as a conditional probability, rather than the 
probability density associated with a pure state. According to Ford et al. [331I34|. this 
normal distribution is associated with two measurements at two different times. For 
each degree of freedom, the overall time-dependent spreading of the wave packet can 
be written as 

wf{t)=(Tl+a-^{t) + s,{t), (14) 
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with i = x,y. Here the quantum contribution to the spreading is given by 

[x{0),x{t)]^ 



4(70 ' r\ 



and Si{t) is the mean-square displacement (MSD) along each direction, 

S,(t) = ({x(t)-X(0)}2}, 

Sy{t) = {{y{t)-v{0)Y). ^'''> 

In order to understand the effects of the measurements on the quantum system, 
apart from the survival probabihty seen in ([2]), one can also analyze the decay rate of 
its initial state by monitoring the ratio between the probabilities at t and t — 0, and 
evaluated at {x,y) = (0,0) Thus, in the case of the initial Gaussian state ([H]), 
we find 

Wx[t) Wy(t) 

which, taking into account (jl4[) . can be explicitly written as 

2 2 

^^'^ ^ ^ ^ al + s{t) + <y'{ty ^^^^ 
since the spreading along each direction is obviously the same. Accordingly, this 
calculation reduces to the simple evaluation of the quantum spreading ai{t) and the 
MSD Si{t) along each direction (i = x,y). The same time-dependence is obtained if 
Gaussian integrations are carried out in the corresponding ratio instead of evaluating it 
at (x, y) = (0, 0). The quantum spreading depends on the commutator of the position 
operators at two different times, from which 

'^'W = F^V^ 'J^'C^*)' (19) 

with ^{rjt) — 1— exp(— ?7t). For each degree of freedom, one obtains the same spreading 
because the Gaussian initial with is assumed to be the same for both directions. 
Analogously, the MSD along each direction takes the form 

s^{t) = ^ t^H{7^t;T). (20) 
7rm 

Here, the function 

has an analytical solution, which allows us to express (PU)) as 

(22) 



.(0 ^ ''"^ 



t - -<^{'nt) 

V 



mri 

In order to analyze ([T8| now for high surface temperatures, we will consider two 
time regimes: rjt (short time) and Tyi 3> 1 (long time). In the short-time regime, 
ryt <C 1, the quantum spreading goes like t^, according to 

a\t) . ^ t\ (23) 
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which corresponds to the wave packet spreading in the absence of dissipation. 
Analogously, the MSD also goes like t^ according to 

s{t) « M: t\ (24) 
m 

where the prefactor is the thermal velocity in two dimensions. Thus, the overall time- 
dependent spreading can be expressed in a more compact form as 

w\t)^al + {v^)t\ (25) 

with 

In this short-time regime (friction- free motion or ballistic regime), we have usually 
that {v'^)t^ < al (the wave packet has not spread too much compared to its initial 
spreading which can be assumed to be arbitrary large) and therefore becomes 



Rit) « 1 - ^t^ + e^("'>tV'^o' ^ (27) 

which is the standard short-time, t^-behavior usually associated with QZE. Notice 
that the effective decay rate [rjesir) — t\J {v'^) /(Jq) is independent of the total friction 
ry and, in particular, of A. Hence the dynamical system displays a friction- free motion 
with a quadratic time behavior of the decay rate. By replacing t by t/N (with 
N oo), some authors [35] show that the factorization given in ^ also applies for 
R{t) and claim that in the absence of friction, the QZE always holds. However, after 
Kofman and Kurizki, for Ohmic friction or linear spectral density, the QZE is not 
expected to occur (see Section [5]). In the indirect measurement scheme used here, the 
dynamical system and the reservoir are entangled at all times except for such a regime, 
since 77 = 7 -I- A is not playing any role yet. Thus, the QZE does not hold because 
any indirect measurement through A has been carried out. Moreover, decoherence is 
absent in this short-time regime since it is a free-motion regime. 

In the long-time (or diffusion) regime, 77^ ^ 1, the quantum spreading in both 
directions is time- independent, with 

^2 

cr^(^^oo) » ^ , ^ (28) 

On the contrary, the corresponding MSD is linear with time, according to 

s{t) « t, 29 

mq 

and Einstein's law is fulfilled, with the diffusion constant given by 

D=^-^ (30) 
mq 

in both directions for an isotropic surface. Thus, the full time-dependent spreading 
can be expressed as 

w\t)^al+ ^ , 2 +^^- (31) 

Assuming ctq ;3> s{t) ^ cr^(t) for a certain time within the long-time regime, R{t) can 
be approximated by 

Rit) « 1 - t. (32) 
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Therefore, by increasing the coverage or cohisional frequency, A, according to ([5^ . 
the decay rate decreases, but obviously one cannot speak about the QZE, because the 
system has aheady relaxed dynamically and decoherence is already manifested. 



3.2. Low and zero surface temperatures 

At low surface temperatures {rj~^ <^ tc ^ ft/?), the noise autocorrelation function is 
complex [36^ and depends on the ratio between the interval of the two times and tc 
(colored noise) . In general, the noise function acts like a driving force and the surface 
dynamics is better described within the generalized Langevin framework as 

x{t) + / dt'ri,{t - t')x{t') = 5F^{t), 

(33) 

m+ dt'Tjy{t~t')y{t)^SFy{t), 



where r]x{t) and riy{t) represent the time-dependent frictions or memory functions 
along each direction. The one-dimensional expression for the noise function can be 
found in the literature [32j . If we assume the surface is isotropic and both frictions are 
Ohmic, ([55)1 reduces to ([TT|). Thus, the corresponding quantum mechanical process 
is not a Markovian process in the customary sense of the term [32] ■ The quantum 
spreading is the same as before, i.e., given by (jl9p . since it is independent of the surface 
temperature. However, the MSD along each direction, given by (PH)) . now reads as 



mr] 

with 



t - -<^{r]t) 
V 



(35) 



27rn 

being the so-called Matsuhara frequencies^ which come from the Taylor series expansion 
of the coth- function in the integral of (|2T|) . At high surface temperatures, p4|l reduces 
to (|22p . The sum over n plays an important role only at very low surface temperatures. 
Nevertheless, at zero surface temperature, the sum disappears, since the coth-function 
becomes unity and H{rit;0) acquires a different dependence on time. In this case 
(T = 0), the corresponding MSD expression reads as 

s,it) ^ ^ (7^ + Inyyt - i [e"^' Ei{~rjt) + e'"^* Eiirjt)] ] , (36) 
nmrj 2 J 

where 7e = 0.577 is Euler's constant, and Ei{—rit) and Eiijjt) are the exponential 
integrals [37]. The environment no longer transfers energy to the adparticle due to 
the zero-point motion of the surface oscillators. 

A similar dynamical analysis can also be carried out in terms of the short- and 
long-time regimes. The quantum spreading in both dimensions remains the same and 
is given by (|23|) and (|28|). respectively. However, the MSD is different due to its 
temperature dependence. Thus, in the short-time regime, T^i <C 1, we find 

sit) ^^-RLf^^e (37) 

m TTTO 

which depends linearly on 77. Analogously, at zero surface temperature, we have 

s(0« — (\-iE-\nTit\, (38) 
7rm V 2 I 
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where we observe essentially the same friction and time dependence, since the presence 
of the friction in the log-function is much weaker. On the other hand, in the long-time 
regime, it can be shown that 

, , 2kBT 2h 1 

s(t)«— , (39) 

mrj nm rj -\- vi 

which, at zero surface temperature, becomes 

« ^(7£; -I- In Tyt). (40) 

Taking into account these behaviors, the following conclusions are drawn at low and 
zero surface temperatures: 



• Short-time regimes. The same formal expression as that given by (pT]) is obtained 
for R{t), except for now (u^) depends linearly on ry, i.e. on A, even at zero 
temperature. Thus, by increasing the coverage or A, R{t) decreases (contrary 
to the manifestation of the QZE) . 

• Long-time regime. The R{t) function can be expressed as 

Kit) « 1 - (41) 

where the MSD goes with A^^ according to ([M]) and (PH)) . Thus, by increasing 
A, the decay rate decreases. However, as before, this manifestation cannot 
be attributed to the QZE, because the relaxation process has already been 
established and decoherence dominates the diffusion process. The AZE is not 
found either. 

For ideal measurements, at zero surface temperatures, Dorlas and O'Connell [55] 
showed that the QZE is characterized by small r]t values, whereas the AZE by large 
values of rjt. According to the previous analysis, though, we conclude again differently 
due to the non-selectiveness or indirectness of the measurements we have considered 
on the time evolution of the adparticle. 



4. The quantum intermediate scattering function 

4-1. Flat surfaces 

For flat or weakly corrugated surfaces where the thermal energy, /cbT, is greater than 
the diffusion barrier, the quantum motion of the adparticles can be exactly described 
within the Heisenberg picture and the scenario of the ISA model [8ll38] . As mentioned 
before, this study can be carried out by means of the so-called intermediate scattering 
function. This function provides us information about the time decay associated with 
the position of the adparticle, which is initially located at R(0) and, after a time i, it 
will be in a position R(i), as it can be inferred from the autocorrelation function ([1]). 
In spin-echo experimental techniques, this function is proportional to the polarization, 
with the real and imaginary parts of the latter being observable magnitudes [UE] ■ 

In general, an exact, direct calculation of /(AK,i) results difficult to carry out 
due to the non-commutativity of the adparticle position operators at different times. 
The product of the two exponential operators in ([1} can be evaluated according 
to a special case of the Baker-Hausdorff theorem {disentangling theorem), namely 
glgB ^ g[l,B]/2gA+B^ ^i^i^ X = -iAK • R(0) and 13 = iAK • R(t), which holds 
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whenever the corresponding commutator is a c-number. When no adsorbate-substrate 
interaction potential is assumed, the formal solution of (|lip can be written as [51I51I55] 

R(i) - R(0) + $(?7t) + — / <^{rit-rit')5¥{t')dt', (42) 

77177 ^"n Jo 

where P(0) is the initial adparticle momentum operator. The commutator of the 
position operators at two different times is a c-number, since the commutator 
[R(0),5F] is zero. This can be justified considering the Caldeira-Leggett model, 
where the two noise functions depend on time, on the initial positions and momenta 
of the harmonic oscillators associated with the two baths, and therefore only on the 
adparticle initial position [SHE]. Thus, /(AK,t) can be expressed as 

/(AK,t) = /,(AK,i)/,(AK,t), (43) 

i.e., as the product of a quantum- mechanical intermediate scattering function, Iq, and 
a classical-like one, Ic- 

Assuming a Markovian scenario, Iq reads as [51155] 



/^(AK,^) = exp 



-^<i>(r,i) 



(44) 



where Er = /i^AK^/2?7t, is the adsorhate recoil energy, which becomes less important 
as the adparticle mass increases. The /^-factor is a time-dependent phase with its 
argument decreasing with the total friction and therefore with the coverage; its time- 
dependence only comes through the ^{rjt) function. At short-times {rjt <C 1 or 
hp <^ ^ivt) ~ Vt s-nd the argument of Iq becomes independent of the total 

friction, thus increasing linearly with time. On the contrary, in the asymptotic time 
limit {rjt ^ 1 or ft/? ^ V~^)i this argument approaches a constant phase. Regarding 
the classical-like factor, 7^, in the Gaussian approximation, it takes the form |5| 

7c(AK,i) ~ f^{t-t')cAt')dt' ^ g-AK^[/(t)+s(t)]^ (45) 

where Cy (t) is the velocity autocorrelation function, which can be expressed in terms 
of the time-dependent functions f{t) and g{t), given by 

fit)-(^-^T;^)i^'"' + Vt-l] (46) 



mj3rf 2mri 



and 



9it) = —l^ 2_,.2, ■ (47) 



2 Vne~''* - rie'""- + rj 

The total intermediate scattering function p3|) can then be expressed as 

/(AK, t) = e-x'l'^'^t-^M]^-AK^git)^ (48) 

with — AK^(t;^)/77^ and a = 1 +ihl3r]/2, {v^) = l/mfi being the thermal square 
velocity. The recoil energy is included in the imaginary part of the product x^oi* i 
which disappears when h ^ Q. Equation P5)) is the generalization of the intermediate 
scattering function for the quantum motion of interacting adsorbates in a flat surface. 
The dependence of this function on AK^ through the so-called shape parameter, x^ 
is the same as in the classical theory [38j. From now on, we will just focus on p8)) . 
which will be analyzed in terms of the QZE and AZE considering the same two time 
regimes as before (short- and long-times). 
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(49) 



where the real part has also the typical t^-behavior characteristic of the QZE, 



which is friction free. As before, not always one has a quadratic time dependence, the 
QZE holds. Indeed, with Ohmic environments, as in our case, the QZE does not take 
place. The imaginary part, in turn, is linear with time, according to 



The short-time regime is related to a Gaussian behavior of the real part of the 
intermediate scattering function, which has been observed in real experiments [2l|3]. 
The time Fourier transform of this function is the so-called dynamic structure factor 
(which is also observable) , leading to Gaussian wings in the corresponding quasi-elastic 
peak. This frictionless motion (no collisions among adsorbates and no interaction with 
the substrate) is the regime where dynamical coherence prevails, i.e., adparticles keep 
their memory of velocity. 

In the long-time (diffusive) regime {-qt 3> 1), p8|) can be written as 



where it is clear that the real part exponentially annihilates any trace of the oscillations 
associated with the phase depending on the recoil energy. Here, the Markovian 
character of the evolution is very noticeable and the regime of coherent dynamics 
is completely lost, which is a typical trait of diffusion regimes. The corresponding 
dynamic structure factor displays a Lorentzian shape around zero energy transfers of 
the quasi-elastic peak. This gradual change of shape in the dynamic structure factor is 
known as the motional narrowing eSect |39fl41) . Actually, strictly speaking, the quasi- 
elastic peak is a mixture of both shapes: Lorentzian-like at very small frequencies 
and Gaussian- like at very large frequencies. According to (|5^ . the real part of 
the intermediate scattering function decays slower with the coverage and therefore 
the dynamic structure factor becomes narrower (Gaussian wings) [30l|42]. At this 
time regime, decoherence is already manifested (or relaxation process). In other 
words, the loss of coherence (or loss of quantum features) begins at times greater 
than the coherence time tc- Quantum motion only plays an important role around 
this time and when t ^ tc, the classical intermediate scattering function dominates, 
since the commutator in Iq plays no role and Ic approaches the classical result. The 
lowest surface temperature reached in a typical ^He spin echo experiment is around 
100 K, which means tc — 0.074 ps. For neutron spin echo measurements, the surface 
temperature can be slowed down up to a few Kelvin, leading to a better analysis for 
quantum coherence. 

Quasi-elastic He-atom scattering measurements |43| constitute an interesting 
example worth mentioning, since they provide a clear evidence for a two-dimensional 
free gas of Xe atoms on Pt(lll) at low coverages {9 = 0.017), low incident atom energy 
{Ei — 10.15 meV) and surface temperature of 105 K along the (100) direction. The 
corrugation of the surface is assumed to be negligible at that surface temperature, 
which manifests in a dynamic structure factor displaying a Gaussian-like line shape. 
Apart from quite different mean free paths extracted from the experiment [43] and 



Re{I{AK,t)}^l---^t\ 



(50) 



Im{/(AK,t)}«-^i. 



(51) 



I{AK,t) « e-x'''*e*■^'-*/^ 



(52) 
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from our own analysis |41) . one would be tempted to consider this system an ideal 
candidate to observe QZE, since the intermediate scattering function is quadratic in 
time for a quite long time interval. However, due to the heavy adsorbate mass, the 
recoil energy is expected to be really small and a classical diffusion motion should 
be good enough to describe this two-dimensional free gas. The classical intermediate 
scattering function for a flat surface is 

/classical (Ai^, t) = CXp [-X^ (e^"* + 7]t - l)] , (53) 

which comes from considering h = in the complex number a and g{t) = in ((48]). 
When r/t ^ 1, the <^-behavior is again obtained 

/classical (Ai^, t) w exp{-6V), (54) 

with S = AK'^fcBT/2m. This behavior with time remains until the adparticle has 
travelled a distance of the order of the mean free path and it starts to lose its coherence. 
But again, for Ohmic friction, the QZE does not hold. 



4-2. Corrugated surfaces 

Surface diffusion of adsorbates on corrugated surfaces can be very often described in 
the Markovian approximation and ISA model by the standard Langevin equation f8) 

mR= -m77R+ Vy(R) (55) 

where T^(R) is the adiabatic interaction potential between the adsorbate and the 
surface. Again the formal solution of (1551) is given by 

R(t) = R(0) + $(7yt) + — / ^i'nt-7^t')[VV(R{t')) + 6F{t')]dt'. (56) 
mrj mrj Jq 

The presence of the adiabatic force introduces an additional commutator, 
[Ro, Vy(R(t))] = {ih)dVV{'R{t))/dPo, where the dependence of this force on the 
initial state (R(0),P(0)) is through R(t). Within a quantum Markovian framework 
[38] . this commutator is very small or negligible. Thus, in the quantum intermediate 
scattering function the Iq factor is the same as in the case of a flat surface, 

given by (|44p . determining again the very short time, friction free behavior of the 
real part of the total intermediate scattering function. The classical- like factor, Ic, 
deserves special attention, since exact quantum Langevin calculations for corrugated 
surfaces are in general, prohibitive. Only the long time regime can be analyzed in 
simple terms. Quantum corrections at different coverages and temperatures in the 
activated surface diffusion of Na atoms on Cu(OOl) have been analyzed [38]. We have 
shown that within the Gaussian approximation, Eq. (|45l) , and by assuming a velocity 
autocorrelation function given by 

C,{t) = ^ e"''* cos((I;i + S), (57) 



Ic becomes 

Ic{^K,t) 



g-xr/(",t) 

^ -> \ m-\-7i 



^-^ mini 

m,n— 



{■m-\-n)fjt-\-i{m~n) {Cjt-\'S) 

(58) 
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where 



and 



/(w, t)=Ai+ A2t + iae^^e-*''-*'^)* + Aie-'~^e-'-^+'^^\ (59) 
~ [ZfjCj sin S+{oj^ -fi"^) cos 6) 

^1 = WT^' ' ^ ^ 

A, = v'ivcosS-u^sin~6)^ ^^^^ 
is = ;t7-^^, (62) 

-2 

= (63) 
2(77 + zo;)^ 

Usually, the values of the parameters f^, <I; and (5 are obtained from a fitting of the 
numerical results issued from solving the standard Langevin equation with periodic 
boundary conditions to (1571) . Moreover, Xi(^K) is a generalized shape parameter 
proposed by the authors to be '44' 



where riy(AK) represents the inverse of the correlation time and is expressed as 

r,(AK) = i.^Pj[l-cos(j-AK)], (65) 
j 

with V being the total jump rate out of an adsorption site and Pj the relative 
probability that a jump with a displacement vector j occurs. 

Notice the linear dependence on time in / because of the parameters f), uj and 
5 are time-independent. This leads to an effective decay rate in ([55]) given by x^^ii 
which accounts for the diffusion and causes the intermediate scattering function to 
vanish at asymptotic times. This scattering function also provides information on the 
low vibrational motion or T-mode. This fact is better appreciated in the dynamic 
structure factor which is the time Fourier transform of the intermediate scattering 
function, 



-X?Ai °° ^_^\m+n 

-T ^-^ m! n\ 



xf^2 + (m + n)f] 



—■ (66) 

[a; — (to — n)CjY + [x^^2 + (™ + n)f]Y 

This general expression clearly shows that both motions (diffusion and T-mode) cannot 
be separated at all. The quasielastic or Q-peak is formed by contributions where 
TO = n. Analogously, the contributions to the T-mode or T-peaks come from the 
sums with n ^ m. If the Gaussian approximation is good enough, the value of r) 
will not be too different from the nominal value of 77 and, therefore, both peaks will 
display broadening as rj (or A) increases. As an illustration, in Fig. [T]the widths of 
the Q-peak for the diffusion of Na atoms on Cu(OOl) at two different coverages and 
two different surface temperatures are plotted. As clearly seen, the experimental [JS] 
and theoretical [301112] values of the Q-peak width display broadening. Taking into 
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Figure 1. Diffusion of Na atoms on Cu(OOl). Numerical (triangles) and 
experimental (circles) dependence of the broadening of the Q-peak on AK , at 
200 K (open circles and triangles) and 300 K (full circles and triangles) along the 
azimuth [100] for two different values of the surface coverage: (a) B = 0.064 and 
(b) e = 0.106. 



account the preceding discussion, this imphes a faster decay rate of the intermediate 
scattering function with the coverage which should be attributed to a manifestation 
of the AZE. 

5. Conclusions 

The characteristic quantum short-time behavior of the observable magnitude (t^-law) 
is always present and, contrary to what is widely accepted, this behavior is not always a 
clear-cut characteristic of the QZE. In general, it depends on the type of measurements 
carried out during the time evolution of the system. This is true, at least, within this 
context (activated surface diffusion of interacting adsorbates) , where the evolution of 
the intermediate scattering function in this time regime is governed by friction-free 
(ballistic) motion. Moreover, within a Markovian framework, the AZE in surface 
diffusion by corrugated surfaces has been reported by experimental and theoretical 
works, but it has not been recognized as such. 

In this work, we have left on purpose for a future investigation several important 
points in the surface diffusion and how they are related to the QZE and AZE. First, 
the diffusion mediated by tunneling. Second, the non-Markovian character of the 
environment and therefore the role played by the sub-Ohmic and supra-Ohmic regimes. 
Third, at very short times, the adparticle is mainly inside any of the potential wells 
of the corrugated surface and the corresponding quantum motion is mainly intra- 
well within a harmonic potential. This simple model allows us to understand the 
physics associated with the so-called T-mode (or frustrated translational mode). 
Recently, QZE and AZE have been shown to occur for the undamped harmonic 
oscillator in terms of the survival probability [28]. A similar analysis could be 
carried out in terms of the intermediate scattering function. Fourth, and finally, a 
true quantum mechanical calculation of the surface diffusion by using the so-called 
stochastic Schrodinger equation of the quantum theory of open systems |llj could 
also be carried out. Work in these directions is currently in progress. 
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